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ABSTRACT 

We show that the excellent optical and gamma-ray data available for 
GRB 080319B rule out the internal shock model for the prompt emission. The 
data instead point to a model in which the observed radiation was produced 
close to the deceleration radius (~ 10^^ cm) by a turbulent source with random 
Lorentz factors ~ 10 in the comoving frame. The optical radiation was produced 
by synchrotron emission from relativistic electrons, and the gamma-rays by in- 
verse Compton scattering of the synchrotron photons. The gamma-ray emission 
originated both in eddies and in an inter-eddy medium, whereas the optical radi- 
ation was mostly from the latter. Therefore, the gamma-ray emission was highly 
variable whereas the optical was much less variable. The model explains all the 
observed features in the prompt optical and gamma-ray data of GRB 080319B. 
We are unable to determine with confidence whether the energy of the explosion 
was carried outward primarily by particles (kinetic energy) or magnetic fields. 
Consequently, we cannot tell whether the turbulent medium was located in the 
reverse shock (we can rule out the forward shock) or in a Poynting-dominated 
jet. 

Subject headings: radiation mechanisms: non-thermal — relativistic turbulence 
— gamma-rays: bursts 



1. Introduction 

Two major unsolved questions in the field of gamma-ray bursts are: (i) the mechanism 
by which the energy in relativistic jets is converted to random particle kinetic energy, and 
(ii) the radiation process by which the particle energy is converted to gamma-ray photons. 

A number of ideas have been proposed for the conversion of jet energy to particle 
energy (cf. Piran 1999, 2005; Meszaros 2002; Thompson 1994, 2006; Lyutikov & Blandford 
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2003; Zhang 2007). Most popular among these is the so-called internal shock model (Piran, 
Shemi & Narayan 1993; Rees & Meszaros 1994; Katz 1994), in which different parts of the 
relativistic GRB jet travel at different speeds. Faster segments collide with slower segments 
in shocks, and a fraction of the jet kinetic energy is converted to thermal energy. Gamma- 
rays are then produced by cither the synchrotron process or the synchrotron-self-Compton 
(SSC) process (e.g., Piran 2005, Meszaros 2002, Zhang 2007). Another model is the external 
shock model (Dermer 1999) in which gamma-rays are produced via the synchrotron process 
in an external shock driven into the circumstellar medium by the GRB jet. Difficulties with 
this model have been pointed out by a number of authors (cf. Piran 1999). 

The excellent data obtained by the Swift and Konus satellites for GRB 0803 19B - 
dubbed the "the naked eye burst" - has provided a new opportunity to investigate the 
viability of these models and to understand the fundamental nature of GRBs. We summarize 
here the main observational properties of this burst (details may be found in Racusin et al. 
2008). 

GRB 080319B lasted for about 50 s and had a burst fluence in the 20 keV - 7 MeV 
band of 5.7 ± 0.1 x 10~^ erg cm~^, which corresponds to an isotropic energy release of 

= 1.3 X lO^'' erg (Golenetskii et al. 2008) for a redshift z = 0.937 (Vreeswijk et al 
2008; Cucchiara & Fox 2008). The time-averaged gamma-ray spectrum during the burst 
had a peak at around 650 kcV. The maximum flux at this energy was ~ 7 mJy, and the 
time averaged flux was ~ 3 mJy. The time-averaged gamma-ray spectrum was measured by 
Konus- Wind (Racusin ct al. 2008) to be F^, oc z/0-i8±o.oi photon energies below 650 keV 
and oc i/-2-87±o.44 higher energies. (The spectrum evolved during the course of the 
burst, as we discuss in §2, and this provides additional information on the radiation process.) 
In the optical band, at photon energies around 2 eV, the peak flux of GRB 080319B was 
V — 5.4 mag or 20 Jy, and the time-averaged flux was about 10 Jy (Karpov et al 2008). The 
optical lightcurve varied on a time scale of about 5 s, while the 7-ray flux varied on time 
scales of ~ 0.5 s. 

GRB 080319B has seriously challenged our understanding of gamma-ray bursts. The 
problems posed by this burst are described in Kumar & Panaitescu (2008) and discussed 
in further detail in the present paper. We show that the observations cannot be explained 
with any of the standard versions of the internal shock model. We find, however, that a 
consistent model is possible if we give up the idea of internal shocks and instead postulate 
that the gamma-ray source is relativistically turbulent (see Narayan & Kumar 2008). 

We begin in §2 by arguing that the radiation in GRB 0803 19B must have been produced 
by the SSC mechanism. Following this, we derive in §3 the basic equations describing a GRB 
that radiates via SSC. In §4, we combine these equations with the internal shock model and 
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attempt to explain the data on GRB 080319B. We find that no consistent model is possible. 
In §5, we consider a relativistically turbulent model, again with SSC radiation, and show 
that in this case it is possible to obtain a consistent model of GRB 080319B. We summarize 
the main conclusions in §6. The Appendix discusses the effects of source inhomogeneity 
on our calculation of the synchrotron self-absorption frequency, and the synchrotron and 
inverse- Compton fluxes. 



2. Why synchrotron-self-Compton model? 

Racusin et al. (2008; Fig. 3) and Wozniak et al. (2008; Fig. 4) show that the 7-ray 
and optical lightcurves (LCs) of GRB 080319B have a similar general shape (although the 
7-ray LC is a lot more variable than the optical LC). This suggests that the optical and 
7-ray radiation were produced by the same source. An independent theoretical argument 
in support of this conclusion is given in §4.1. The radiation mechanisms in the two bands 
must, however, be different since the optical flux is larger by a factor ~ 10^ than the 7-ray 
flux extrapolated to the optical band. 

The average spectral properties of GRB 080319B were summarized in §1. Racusin et al. 
(2008; supplementary material - Table 1) also reported spectral flts corresponding to three 
independent time segments of the burst: —2 s to 8 s, 12 s to 22 s, 26 s to 36 s (all times 
measured with respect to the nominal start time of the burst). The peak of the gamma-ray 
spectrum evolved from about 750 keV early in the burst to about 550 keV at late times, 
giving a mean peak energy of 650 keV as mentioned earlier. More interestingly, the spectral 
slope at energies below the peak evolved with time: oc //0-50±0-04^ jy0.i7±o.02^ j^o.ioio.os^ 
during the three time segments. The unusually hard spectrum during the flrst time segment 
unambiguously points to an SSC origin for the gamma-ray emission, as we now argue. 

The hardest spectrum possible with optically thin synchrotron emission is F^j oc u^/'^ . 
The only way to obtain a harder spectrum is to invoke self-absorption, in which case the 
spectrum will switch to F^ oc v"^ below the self-absorption break. However, in order to obtain 
a synchrotron spectrum with a mean spectral index of 0.5 in the band between 20 keV and 
650 keV, we would need to have the self-absorption break at an energy ~ 50 — 100 keV. This 
has two serious problems. First, it is virtually impossible to push the self-absorption break 
to such a large energy with any reasonable parameters for the radiating mediuno. Second, a 



^For the synchrotron self-absorption frequency to be ~ 50 keV and the flux at 650 keV to be ^ 10 mJy, 
the distance of the source from the center of the explosion must be less than 10^ cm. At such a small 
radius the medium would be extremely opaque to Thomson scattering and 7 + 7 — > e+ + e^. Therefore, the 
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spectral break in which the slope changes from 1/3 to 2 would almost certainly be detectable 
in the data and would not be consistent with a single power-law with a slope of 0.5. 

A Comptonization model gets around these difficulties. If the gamma-ray emission 
is produced by Compton-scattering, then any break in the spectrum is not intrinsic to the 
gamma-rays but merely a reflection of a break in the spectrum of the underlying soft photons. 
If the soft photons are in the optical-infrared band and are produced by synchrotron emission, 
then the self- absorption break must be around 0.1 eV, which is perfectly compatible with 
reasonable model parameters. In addition, although the synchrotron spectrum below the 
break would be very hard, viz., Fi, oc z/^, the corresponding segment in the up-scattered 
inverse- Compton radiation would be softer: F^, oc u (Rybicki & Lightman, 1979). Thus, the 
gamma-ray spectrum would break from a slope of 1/3 to 1, which is consistent with the 
observations, especially when we allow for a smooth rollover from one spectral slope to the 
other over a range of energies. 

Why did the gamma-ray spectral slope below 650 keV switch to ~ 0.1—0.2 at later times? 
The likely explanation is that the self-absorption frequency of the synchrotron emission 
dropped to yet lower energies (in the infrared), and so the break in the gamma-ray band was 
pushed closer to, or even below, 20 keV. 

In this discussion, we have assumed that the soft radiation is produced by the syn- 
chrotron process. The alternative is thermal radiation, but this can be ruled out as it 
requires a Lorentz factor of ~ 10^ {6t)~^T^^^'^ to explain the observed flux of 10 Jy@; here 
St is the observed variability time (in seconds) of the optical lightcurve and T5 = T/IO^K is 
the temperature of the source. The parameter T5 cannot be much larger than unity since 
the total energy release would become excessive (> 10^^ erg). A Lorentz factor of 10^ is not 
reasonable either. Therefore, a thermal model for the optical emission is ruled out. 

We thus conclude that the optical photons in GRB 08030 19B were produced by the 
synchrotron mechanism, and the gamma-ray photons were produced by the same relativistic 
electrons by inverse- Compton scattering the synchrotron photons. That is, all the observed 
radiation in GRB 080319B was the result of the SSC process (§4 gives a more detailed 
discussion) . 



emergent radiation would be thermal and no photons with energy > 1 MeV would be able to escape from 
the source. However, Konus-Wind detected ~ 10 MeV photons from GRB 080319B. 

^For a relativistically moving thermal source with a temperature T (in the observer frame) and ra- 
dius i?,, the observed flux in the optical band at a frequency I'op ^ 5 x 10^^ Hz is: fop ~ 2kT{l + 
zYvlj,R^ /{dlc^V^) w 2fcr(l + zYh'lj,{6tfV'^ /d\. Therefore, to explain the observed optical flux of 10 Jy we 

require T - lQ^{5t)-^T;r^''^ . 
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3. Synchrotron and SSC processes for a relativistic transient source: basic 

equations 

We first determine what properties the source of optical emission in GRB 080319B must 
have, assuming that the radiation is produced by synchrotron emission (§3.1). We then 
consider the gamma-ray data and describe the additional constraints they provide (§3.2). 
Detailed application to GRB 080319B is discussed in §4 and §5. 



3.1. Modeling the prompt optical data 

The properties of a synchrotron source can be described by five parameters: B, Ne, F, 7, 
and Te, which are the magnetic field strength, the total number of electrons, the bulk Lorentz 
factor of the source with respect to the GRB host galaxy, the typical electron Lorentz factor 
in the comoving frame of the source (the electron distribution is drie/d'y oc 7"^ for 7 > 7j), 
and the optical depth of the source to Thomson scattering. If the source is at redshift z, the 
peak frequency z/j of the synchrotron spectrum and the observed flux fi at the peak, both 
as seen by the observer, are given by (Rybicki & Lightman, 1979) 

where di is the luminosity distanc^to the source, and and (j)f are dimensionless constants 
that depend on the electron energy distribution index p; for p = 5 (as suggested by the high 
energv spectral index for GRB 080319B), (p^, = 0.5 and (pf = 0.7 (cf. Wijers & Galama, 
199!cl). In all subsequent equations we explicitly use these values for the 0's with the 
exception that for the most important quantities - such as the total energy and the inverse 
Compton flux - we show the dependence of the results on the 0's to indicate how uncertainties 
in Ui and fi affect the final result. 

Throughout the paper we measure frequencies in eV and fluxes in Jy. All other quantities 
are in cgs units, but we use the short hand notation Xn = ^/(lO^cgs) to scale numerical 



•^Technically, the electron index p is a sixth parameter, but since it can be estimated directly from the 
high energy slope of the gamma-ray spectrum we do not count it 

"^The factor (1+z) in the numerator of the expression for flux, eq. 2, is due to the fact that the luminosity 
distance refers to the bolomctric flux whereas we are considering the flux per unit frequency. 

^(fii, = 1.5xp and (f>f = (j)p in the notation of Wijers & Galama. 
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values, e.g., R15 = R/ (10^^ cm) is the radius of the radiating shell with respect to the center 
of the explosion in units of 10^^ cm. Using this convention, the source optical depth Xg and 
the duration of a pulse in the lightcurve 6t (in seconds) are given by: 

re = ^ = 0.5^ or N^,, = 2t^RI, (3) 

6t = ^^^^ = (1-7 X ioS)/?i5r-^(i + z). (4) 

Equivalently, 

r = 130(5t)-i/'(l + z)^/2i?}/'. (5) 
We can combine equations (1) & (2) to eliminate B: 



h = ^ \°'^^ .iV,55(l + zYdll, (6) 



and use equation (3) to replace by Xg. The resultant equation is: 



f. = ^-^''l^ ^'' ^.YRUl + zfdll,^ (7) 

or 

7. = 8li^"^/VVV^i^lf(l + .)V2rf,i/^ (8) 

Here we have defined the quantity 

y ^ 7^e, (9) 

which is closely related to, but smaller by a factor of 2-3 compared with, the Compton- 
Y parameter for GRB 080319B (the difference depends on electron distribution function). 
Substituting eqs. (5) & ([8]) back into eq. (1), we find 

B = {2^hG){5tf''ul"fl"Y-"'R-^'\l + z)-'/'d,,,, (10) 

while equations ([3]) & ([8]) yield 

iV,55 = 2.9 X lQ-'v;"^fl'^Y"^R^,{l + zrUL2s. (11) 



Equations ([3]), ([5]), ([8]), ffTOl) and (|TT1) are solutions for the five basic physical parameters 
of the synchrotron source in terms of the observed variability time 5t and four unknown 
quantities: the source radius i?i5, the Compton parameter F, the peak energy of the syn- 
chrotron spectrum z/j and the synchrotron flux at the peak /j. The last two quantities are 
not independent — they are constrained by the observed optical flux fop (= 10 Jy in the 
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case of GRB 0803 19B) at 2 eV. Depending on whether the synchrotron peak frequency z/j is 
below or above 2 eV, we obtain the following constraint: 

(h{2/ui)y\ V, > 2, 

fop = { (12) 

[fM/2y/^, V, < 2, 

where we have used standard results for the spectral slope below and above the synchrotron 
peak, assuming for the latter that the cooling frequency is close to Ui, as suggested by the 
spectrum of GRB 0803 19B§; for other bursts, where the cooling break is substantially above 
the synchrotron peak, the equations in this paper can be easily modified by replacing p with 
(p-l). 

Another constraint is provided by the synchrotron self-absorption frequency Ua, which 
can be estimated by equating the synchrotron flux at Ua to the blackbody flux in the Rayleigh- 
Jeans limit. Assuming that the electrons that dominate at Ua have Lorentz factor 7^, we 
obtain 

-^meC7. = /(.J = ^(^^^(^^J , (13) 

where primes refer to quantities in the source comoving frame. To convert to the observer 
frame, we use 

u,^u'J/{l + z). (14) 
Combining the above two equations we find 



3_8;3/5^-l/5r3/5^-3/5^-6/5^^ + zy^/^dtl (15) 



As before, all frequencies are in eV and fluxes in Jy. Substituting for F & 7^ from equations 
dS]) and dHD results in 

= 5.1 /f/V-'/'°(5t)-='/i°i?-^/V-=^/20(l + z)-'/'d%l. (16) 

As we discuss in §3.2, we can estimate from the gamma-ray data the ratio rj of the synchrotron 
frequency to the self- absorption frequency Ua- 

V = ^i/^a- (17) 



^The spectrum of GRB 080319B peaked at 650 keV and contained no other break between 20keV and 7 
MeV. Since the coohng frequency cannot be smaller than 650 keV since the spectrum varies as v'^'^ below 
the peak, it must be either larger than 7MeV or close to 650 keV. The former possibility is unlikely since 
it corresponds to a radiatively inefficient system and would increase the energy requirement of an already 
extreme burst. 
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Therefore, this is a third observational constraint on the source properties. 

To summarize, the source of optical emission is described by means of five parameters. 
The pulse duration 6t, the optical flux fop and the frequency ratio r] (eq. [T7]) give three 
constraints. The solution space is thus reduced to a two-dimensional surface. Additional 
constraints are obtained from the gamma-ray data, as we discuss next. 

3.2. Gamma-ray emission via the inverse-Compton process 

We assume that the gamma-rays are produced by inverse Compton (IC) scattering of 
synchrotron photons. The peak frequency of the IC spectrum and the flux fic at the peak 
are related to z/j and fi as follows, 

'y^c ^ fic ~ 3re/, = 'iYlff^, (18) 

where a multiplicative factor of 3 in the expression for fic takes into account the ratio of 
solid-angle integrated specific intensity inside the source and the flux just outside the shell 
(in the source comoving frame). Substituting for 7j from equation ([H]), 

Z/,e6 = 1.9 X 10-2/r'^'^f^'>^'/'^15(l + Z)dl\^, (19) 
/.e-3 = 0.47/f/'z.rVVl/2i?-^l(l + ^)-lrf^28, (20) 

where i/^ce = '^ic/(10^ eV) and /^c-a = /ic/(10"^ Jy). 

Using equation fITI?]) for z/j^ we determine the synchrotron peak flux, 

/, = 3.8 X \^-'v-;^,vlR\^{\ + zfd-jl^, (21) 
and substituting this into equation (fT6|) we obtain Vi. 

V, = 1.2 X 10'r]~''/'u^i,\6ty/'R;,'^'V/%l + zf'\ (22) 
Eliminating z/j from equation (12T]) by using eq. (!22|) . we flnd 

/, = 6.6 X lQ^r,-''"^i^iJtR^^Y-\l + zfdll^. (23) 

Substituting for Vi and fi from equations (122|) & (123|) into equation (fT2l) for the optical 
flux we obtain 

166r?-wz.,^«e(5t)ii?,|y-i(l + z)^dZL ^^>2, 

922 X >i.Tr'-^.Tm)'¥RfY-'¥il + .j'^^rf- < 2. 

(24) 
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We use the optical flux to eliminate one more unknown variable, Y, 



18 



7Ur^-^^fo:/u/J6t)^R,l (1 + z)^428 , > 2, 

Y = (25) 

3 5{p+6) 3_ 5p+18 6-p ^^jg _ _6_ 

[471 X 1.4Ff5 7]~^(^ fop"^' {5t)^Rll''' (1 + z)^d^^^'' , z/^ < 2. 

We are now in a position to express all quantities in terms of four observables, viz., the pulse 
duration 5t, the optical flux /op, the dimensionless self-absorption frequency ?7, and the peak 
frequency of the gamma-ray spectrum z/jc, plus one unknown parameter Ri^. 

The peak frequency of the synchrotron spectrum is obtained from equations fl22|) & fl25|) : 

\.^v^flyum-'^R'J{i+zr'^d%,, V, > 2, 

"'^^ 2 1 3 4 ^^^^ 

2 X lA-^^T^^fi'u^^f^'idtr^R^.^'il + zy^^d]^,, V, < 2. 



Let us define the transition radius Rtr as that value of R for which Ui = 2 eV. From equation 
(I2nD we find 

10 2 _i 1 si 

i?t.,i5 = 0.8r]^fIpi^;j{5t)-Hl + zyUl^,. (27) 
Note that, for R > Rtr, i^i < 2 eV, and vice versa. 

The flux at the peak of the synchrotron spectrum is obtained from eqs. fl23|) & fl25|) : 

f 10 5 -1 1-3 si 

\0.92r)-ifJpi^,J{St)-rR,,r^l + ,)-rdl,„ R < Rtr, 

h={ , 3 . 3. ^^^^ 

and the peak gamma-ray flux is obtained from equations fEUl) . fl2^ . fl2Bl) & fl25]) : 

12.3 r/^2T/o>Jg(5t)vi?^/(l + z^d^l, R < Rtr, 



jc— 3 



2p+4 10(p+2) ^ 2p+5 j_ 3p+4 _ 2E±2 

8.4 X 1.4— r]-^(^f^+' vl^^ {bt)—^Rlt^ (1 + z)l^ d^^^^ , R > Rtr. 



(29) 

Using equations ([23]), ([2nD & (EHD to substitute for Y, Ui & fi, equations (fTUl) fc 
TT]) give 

|555r/-v/,74g(5t)^i?M(l + ^)^c^ZJ8, ^ < 

7^ = <^ (30) 

P+4 



527r7" /o/+3iy//g+'^ (5t)^i?i'f (1 + z)^cii2T^'^ . i? > i?, 



tr. 
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4.5G V-foJ, ^^6^ (5t)--i?r5"(l + ^)"^^/28, R < Rtr, 

B = ^-U (31) 



P+5 „_-, P+15 



17-p 



[e.l X 10-3r7-i/o7^|(<5t)fi?S(l + ^)^t^Zl8, ^<^. 



AT, 



e55 



5(p+4) 1— p+4 3p+12 g 2_ 

5.9 X 10-3 r/" zy.^2gP+« (5t)55+6i?if (1 + z)^d^^+' , i? > Rtr. 



(32) 

The energy in the magnetic field as measured in the GRB host galaxy rest frame is 
Eb = B'^R^/2. Using equation (13T|) . this can be shown to be 



1.2 X lO^^erg ^^fj^ m~^R;,^ (1 + zr^dl,,^-^ R < R,r, 



E 



B 



2p+10 „_-, 2p-6 



p-1 ^P-" 17-p 



(33) 

The energy in the charged particles (e^) that produce the optical and gamma-ray emission 
is E^ = Ne-fiTnieC^. From equations dS}, ([30]) & (I32D this is 



1 4 X lO^^erg iy,^{St)i R^,{1 + z)"^ dJ,{ct>J cl>f)h ^ R < 

I (5P+25) 2_ 2£±Z 1 2£+9 p+8 ^ 1 

[3 X lO^^erg v'^^U^'Kt" i^t)^ Rw' (1 + ■2)^c?l28^' ^ ^ > Rtr- 

(34) 

We should also consider the energy in the protons arising from the bulk relativistic motion 
of the shell. However, to estimate this quantity we need to make some assumption regarding 
the composition of the fiuid in the shell, whether it is primarily an e~^e~ plasma or a p'^e~ 
plasma. In the former case the energy in protons is negligible, while in the latter case the 
energy is Ep = NeViripC^ ~ few x E^ (taking 7^ of order a few hundred). 

In the observer frame, the synchrotron cooling time tgyn is 

tsyn = (7.7 X lO's) (1 + z)T~'B-^jr\ (35) 
which, using equations (15U1) and ([SI]), can be written as 

(620s) r]-'f-„^ul,{6t)Rl,{l + z)' dll,, R < Rtr, 



E,, 



,5/2 y^^^'^j '/ J op 



T- 3P+16 4-p £+24 :,9.-r, li_ 



/ Ij295s) r7-OTy/„/+V^2p+6 ^ ^^_^^,+3^ ^ ^ 

(36) 
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In addition to the loss of energy via synchrotron radiation, electrons also lose energy through 
IC scattering of the local radiation field. We calculate the photon energy density in the source 
rest frame from the observed bolometric luminosity Lots, and use this to estimate the IC 
cooling time tic in the observer frame: 

Strictly speaking we do not know the true bolometric luminosity, so the above estimate of 
tic is an upper limit to the actual inverse-Compton cooling time. Using equation fl30|) we can 
rewrite tic in the following more useful form: 



Uc = ^\ ■ (38) 

All the results obtained so far are general and could be applied to any GRB that has 
the required data. We now consider the implications for the naked eye burst GRB 080319B. 



4. Application to GRB 080319B: Ruling out the internal shock model 

The relevant observational parameters for GRB 080319B are: z = 0.94, (1^28 = 1-9, 
fop = 10 Jy, z/jc6 = 0.665, p = 5, and 5t ~ 1 s (from the gamma-ray variabihty). Moreover, 
we estimate thatlll the time-averaged -q = Vi/va for this burst was about 25 (because fv oc 
j^o.iSito.oi between 20keV and 650 keV), whereas the initial value of rj was ~ 10 as oc //0-5±o.o4 
during the first 8 s. Scaling all quantities to these values, the transition radius Rtr becomes 

10 2 _1 ^ 

i?,,,i6 = 12VJ4iZ/-|8(5t)-3. (39) 

If RiQ > Rtr,m, then Ui < 2 eV and the optical band is in the steep decaying part of the 
synchrotron spectrum, above the synchrotron peak. If Riq < Rtr,i6, then the synchrotron 



^The spectral indices are obtained by fitting the data with the Band function, which gives the asymptotic 
value for the low energy and high energy index. It should be noted that the IC spectrum below ^ "^^ifva 
is /,y cx V and between this frequency and the peak of vfi, at i^fvi the spectrum changes from f^(xvio 
v^^ . Therefore, somewhere in between these two frequencies the index would be 0.2. The Konus data for 
GRB 080319B found the spectral peak to be at 665 keV and the observations extended down to a minimum 
photon energy of 20 keV. Therefore, the IC spectral index of 0.2 at 20 keV requires Vi/va ^ 25. 



- 12 - 



peak is above 2 eV, and the optical band is in the F^, oc z/^/^ part of the synchrotron spectrum. 
Since the prompt optical emission in GRB 080319B was exceptionally bright, it is likely that 
the peak of the synchrotron spectrum was fairly close to the optical band. This suggests that 
RiQ must be within a factor of a few of Rtr,i6- According to equation (jSHD, Rtr ~ 10^'' cm, 
which is orders of magnitude larger than the radius at which internal shocks are expected. 
In fact, it is comparable to the deceleration radius of the jet. 

From the results described in §3, we obtain the following numerical results for the 
relevant parameters in GRB 080319B: 

r = 572{6t)-^^^R\i\ (40) 

(77 vU Qi 4.8 m ^^Rlh R< Rtr , 

[251 vitf'^ii^ksm^^Rfl, R > Rtr, 

[(1.1 X 10^1 vilLpl'^lASt)^Ri, R < Rtr, 

Ne = { (42) 

1(3.4 X 10^1 Qi^lkm-^Rll R > Rtr, 



80 9 22 o 10 

~3 , 



2.5 X 10-3r^,.f /op,>J5.8(^i)^^16, R < Rtr, 

Y = I (43) 

1 _55 _ 3 43 11 J_ 

[7.5 X lo-'vi.t&lsm^^Rll R > Rtr. 



40 24 _26 . _22 

(2 X lO^^erg) Vi^AUA^A^tyrR;^^ , R < 

Eb= { ^ ^ ^ (44) 

(2.3 X lO^^erg) vLfopAIsi^tf^RL R > Rtr, 

_65 _6 17 3 23 

' (3.4 X lO^^erg) r^-f Ul^^^ksmi R,l , R < R,,, 

Ee = { (45) 

_50 _1 17 1 15 

,(3.3 X lO^Oerg) ViT Lp^i^^^ksi^tys R,%, R > R.^, 



(7.4 X 10-^s) vLlfapA^LA^ml,, R<R 



■tr , 



tsyn = < (46) 

_35 _7 31 -I 29 

f2s) virUi'yi:,.sm-eRii R > 



2.1s 

'jo6s,52 
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Figure 1 shows the dependences of a number of quantities as functions of the only free 
parameter in the modeh Riq. The lower solid line in the top left panel shows the predicted 
gamma-ray flux /ic-s (based on eq. HHl) and immediately indicates a major problem. If, as 
we suggested earlier, Riq ~ Rtr,i6, then the peak IC flux fi^ that the model predicts falls 
short of the observed flux by nearly a factor of 100 

In Appendix A we discuss possible sources of error in our estimate of the IC flux. We 
show that the uncertainty in /jc, even after allowing for inhomogeneities in the source, is no 
larger than a factor of order unity. The largest error is that we have overestimated Ua by 
a factor ~ 1.5 by not including the expansion of the source during the time it takes for a 
photon to cross the shell (see Appendix A). The effect of this is that fic is underestimated 
by a factor ~ 2.5 due to its dependence on Ua via 77. Even after correcting for this (upper 
solid line in upper left panel in Fig. 1), the theoretically calculated gamma-ray flux is still 
smaller than the observed value by a factor ~ 30. This discrepancy is much too large to 
be overcome by minor adjustments to the model. We thus conclude that the internal shock 
model with Riq ~ Rtr,i6 is ruled out for GRB 080319B. 

One way to mitigate this problem is to select values of Riq that are either very much 
smaller or very much larger than Rtr,w- However, as Fig. 1 shows, we need to modify Riq 
by a huge factor, which immediately leads to other problems. 

The top right panel indicates one of the problems we face. This panel shows the total 
isotropic energy of the source {Eb + Ee) in units of 10^^ erg. We see that shifting Riq 
substantially away from Rtr,i6 causes the total energy to become unphysically large. A 
reasonable upper limit to the total energy is 10^^ ergj, which corresponds to an energy of 
~ 2 X 10^3 erg in each spike in the gamma-ray lightcurve The energy estimates in equations 
(jH]) and fHSl) refer to the latter and the limit is shown by the dotted line in the top right 
panel. We see that the energy constraint restricts Riq to lie within the range 4.4 — 14.6. 
Within this range of Riq, we have approximate equipartition between Eb and Ee (see the 
dashed line), which is desirable, whereas choosing other values of Riq would cause large 
deviations from equipartition. 



^The isotropic energy release for GRB 080319B in the 20keV — 7MeV band was 1.3x10^^ erg. The 
radiative efhciency for GRBs varies widely from burst to burst but is generally larger than ~ 10% (Panaitescu 
& Kumar, 2002). Therefore, the total energy release in GRB 080319B is expected to be < 1.3xl055erg, and 
so we take {Eb + E^) < 10^^ 
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Fig. 1.— Top left: Shows the IC flux in GRB 0803 19B at 650 keV predicted by the internal 
shock model. The lower line is from eq. fl48p and the upper line corresponds to a factor of 2.5 
larger flux to allow for the expansion of the source during a photon crossing time, which was 
not included in calculations presented in §3 and §4 (see the Appendix for details and for a 
discussion of uncertainties in fic and z/q). The shaded band is the region of the source radius 
i?i6 that is favored by various constraints (see text for details). For these values of Riq, the 
predicted IC flux falls short of the observed gamma-ray flux (horizontal dotted line) by a 
factor >30. Top right: The solid line shows the total isotropic energy in a single variability 
spike in the gamma-ray lightcurve. The dotted line shows the maximum energy allowed. 
The dashed line shows Eb/E^. Bottom left: The solid line shows the cooling time in units 
of the variability time (Is) - which should be > 1 since the low energy spectral index was 
positive - and the dashed line shows the Compton Y. Bottom right: The solid line shows 
the bulk Lorentz factor F and the dashed line shows the typical Lorentz factor of electrons 
7i- 
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The bottom left panel in Fig. 1 shows another set of problems. Given the huge luminosity 
of GRB 080319B, we expect the source to be radiatively efficient, which means that tcooi must 
be comparable to the variability time tvar ~ 5t ~ 1 s. We see that the cooling time tcooh 
calculated according to 

1 _ 1 ^ 1 ^^g^ 

^cool ^syn '^ic 

is within a factor of 10 of the variability time only for models with Riq in the narrow range 
2.7 — 6.9. (The variation of tcooi with Riq is extremely steep, so the condition tcooi ~ ^var is 
very restrictive.) Combining this constraint with the one we obtained earlier from the total 
energy, the allowed range of Riq is limited to 4.4 — 6.9, shown as the hatched vertical band 
in the top left panel in Fig. 1. We also see from the bottom left panel that the Compton Y is 
less than 0.1 for most models, and extremely small, Y -C 0.1, for small values of Riq. Since 
Y determines the fraction of the source luminosity that comes out in gamma-rays, and since 
GRB 080319B (indeed, any GRB) is a strong gamma-ray source, it seems unlikely that Y 
could be this small. 

Finally, the bottom right panel in Fig. 1 shows the dependence of the bulk Lorentz 
factor r and the random electron Lorentz factor 7^ on Riq. Models with Riq ~ Rtr,i6 predict 
reasonable values ~ 10^ — 10^ for both Lorentz factors, but models with very different values 
of RiQ predict either unusually low or unusually high values. 

In summary, all the indications suggest that the optical and gamma-ray radiation in 
GRB 080319B were produced at a radius R ~ few x 10^^ cm — 10^^ cm. But at this radius, 
the internal shock model predicts a negligibly small gamma-ray flux. We are thus forced to 
conclude that the internal shock model, at least in its standard form, is definitely ruled out 
for GRB 080319B. 



4.1. Other versions of the internal shock model 

We now consider whether we can get around the above difficulty by modifying the 
internal shock model. We begin by noting that, as long as the gamma-ray emission is IC 
- something that is required by the low energy spectrum fi^ oc u^'^ at early times (§2) - 
and the seed synchrotron photons are produced in the same source as the 7-ray photons, 
equation fH5]l is valid. This equation predicts an unacceptably low flux in the gamma-ray 
band. Therefore, if we wish to save the internal shock model, we must give up the assumption 
that all the radiation came from the same region of the source. 

Let us assume that the seed photons for IC scattering are produced by the same source 
that gave us the optical flash. We will call this the optical region of the source. Let us 
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assume that these seed photons are IC-scattered in a different region, the gamma-ray region. 
We now show that the electron Lorentz factors ji in two regions are very similar. 

Let us suppose that 7j in the optical region differs from that in the gamma-ray region. 
Then, the self-IC radiation from the optical region will introduce a second IC component in 
the observed spectrum, with a peak at a different photon energy. Equation fHSl) is valid for 
any SSC process, so we can use it to estimate the flux in the second peak. If the IC peak 
from the optical region is at a higher photon energy than 650 keV by a factor > 2.5, then 
equation (HHil shows that the flux in this component will be larger than the observed flux 
(note that fic oc ul^^ as per eq. HHl and the observed flux above 650 keV declined as z/~^-^^). 
On the other hand, if the self-IC radiation peaks at an energy much less than 650 keV, the 
magnetic energy in the source will increase very rapidly {Eb oc eq. 111]). Since the 

energy is already close to the maximum limiting value we can accept, this option is also ruled 
out. 

Therefore, the values of 7^ in the optical and gamma-ray regions must be nearly the 
same. This tight relation between the Lorentz factors in the two regions suggests that 
the optical and gamma-ray sources are very likely the same region. Even if they are not, 
the similarity of their parameters means that the large discrepancy in the gamma-ray flux 
discussed previously will survive unchanged. A related idea is that there are two populations 
of electrons with different values of 7^ within the same source. One population is responsible 
for the seed photons and the other for the IC scattering. This possibility can be ruled out 
by the same argument. 

This leads us to consider a model in which part of a shell is magnetized - this is where 
optical photons are produced - and the rest has a much weaker magnetic field (in order to 
avoid overproducing synchrotron flux) but contains about 30 times more electrons in order 
to produce the observed ~ 3mJy 7-ray flux via IC scattering. This situation can arise, for 
instance, when magnetic field decays downstream of a shock front, as suggested in Kumar & 
Panaitescu (2008). However, this proposal suffers from serious problems that these authors 
have pointed out in their paper. The solution requires magnetic field to decay on a length 
scale that is about 5% of the shell thickness or about 10'' plasma skin depth . This scale 
corresponds to no particular physical scale in the system and is quite arbitrary. An even 
more severe problem is that the model cannot account for the shorter time scale variability 
of gamma-rays compared to the optical; in fact, the natural expectation is the opposite in 
this model. 

Note that Fig. 1 indicates an extremely narrow range of R for the radiating medium. It 
is hard to believe that a large number of independent shells ejected from the central source 
would all collide at exactly this radius. In addition, as we noted earlier, the radius R of the 
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source is uncomfortably large for the internal shock model. Both of these features would 
be explained naturally if we assumed that the internal shocks are not between independent 
shells, but rather between successive shells and the outermost shell, which is decelerating 
after colliding with the external medium. This is a variation on the general idea of internal 
shocks (with a strong hint of the forward shock model, see §5.2.2), which at least provides 
an explanation for the radius of the source. However, this model can be ruled out for two 
reasons. As with all other variants, this model cannot explain the magnitude of the IC 
flux unless the magnetic field occupies a small fraction of the shocked shell, about 5% of 
the ejecta width or 10^ plasma skin depths. Furthermore, it predicts that the pulse-width 
should increase with time, which is inconsistent with the observed data for GRB 080319B 
which show, if anything, that the last few pulses in the gamma-ray lightcurve were somewhat 
narrower than the initial few pulses. 

Having considered these and other ideas, we believe that it is impossible to explain the 
observations of GRB 080319B with any reasonable version of the internal shock model. For- 
tunately, there is an alternative model which invokes relativistic turbulence in the radiating 
fluid. We now apply this model to GRB 080319B. 



5. Relativistic turbulence model for GRB 080319B 

The basic kinematic features of the relativistic turbulence model are described in Narayan 
& Kumar (2008). In brief, this model explains the observed variability in GRB lightcurves 
by postulating an inhomogeneous relativistic velocity field in the GRB-producing medium 
(which we refer to as the "shell" because of its shell-like morphology in the host galaxy 
frame). The beaming effect of the turbulent eddies causes large amphtude fluctuations in 
the observed flux. Despite being inhomogeneous, the model is radiatively efficient in the 
sense that the whole medium radiates and the observer receives a fair share of the radi- 
ated luminosity. This important feature, which is a direct consequence of beaming, allows 
the model to overcome the arguments of Sari & Piran (1997) against inhomogeneous GRB 
models. The reader is referred to Narayan & Kumar (2008) for details. 

Since the relativistic turbulence model has a natural explanation for the observed vari- 
ability, equation (jlj) relating the variability time scale 5t to R and F is no longer needed. 
Instead, the quantity R/2V'^c determines the total burst duration t^. We thus have 

t, ~ ^^^^ = (1-7 X 10^ s)R,,T-\l + z). (50) 

Since ~ 50 s for GRB 080319B, whereas 5t ~ t-^ar ~ 1 s, this modification has a rather 
profound effect on the results. 
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In the relativistic turbulent model, we consider turbulent eddies with a typical bulk 
Lorentz factor 7^ in the frame of the shell, and a typical size ~ R/ (7tr) in the comoving 
frame of an eddy. The eddies are volume-filling, so there are ~ 7^ eddies in a causally 
connected region of volume ~ R"^ /T^. We assume that the velocity field of eddies changes 
direction by ~ 2ti on the light crossing time scale ~ R/{c'~ft^). In this case the probability 
that an eddy, some time during its life, will move towards the observer with a velocity vector 
within an angle (7tr)~^ of the hne-of-sight is ~ 7^"^ (Narayan & Kumar, 2008). Therefore, 
over the course of the burst, a given observer will receive emission from 7I eddies, with each 
eddy producing a pulse of radiation lasting a time (see Narayan & Kumar 2008 for details). 

Since GRB 080319 has t^ar ~ t-y/lOO, we infer that 74 ~ 10 for this burst. Note that, at any 
given time, the observer receives radiation from only one eddy on average. 

We assume that the fiuid in the shell consists of eddies and an inter-eddy medium. The 
latter is produced when eddies collide and shock. Let us take the thermal Lorentz factor of 
electrons within an eddy to be '-fu. The thermal Lorentz factor of electrons in the inter-eddy 
medium follows from energy conservation when eddies collide, and is ~ jitji = li- Similarly, 
if we take the magnetic field in the inter-eddy frame to be S, then the comoving magnetic 
field in an eddy is 5/7^^, assuming that the magnetic energy is roughly conserved when 
eddies dissipate. Using these scalings we see that the peak of the synchrotron spectrum (as 
measured in the shell frame) for inter-eddy and eddy emissions are proportional to B^f and 
^lllt'^^'^ 1 respectively. 

Let us take the average number of electrons in an eddy to be A^g^, and the total number 
of electrons in the inter-eddy medium in a volume (i?/r)^ (the volume of a causally connected 

region) to be Ni. For simplicity, let us assume that the total number of electrons in all the 
7^ eddies is of order A^j, i.e., half of the fiuid in the shell is in eddies and the other half is in 
the inter-eddy medium. Thus we have N^^ ~ ^i/lt- 

At any given time, only one eddy will produce beamed radiation towards the observer. 
The peak synchrotron fiux from this eddy is proportional to ~ BN^.^^^^'^V^ ~ BNjT^ /^^^"^ . 
Here we have made use of the fact that, at a fixed observer time, the observer receives radia- 
tion from only a fraction of the electrons in the eddy, ~ Ne/^it- due to the time dependence of 

eddy velocity direction. The peak synchrotron fiux from the inter-eddy medium is ~ BNiT^, 

3/2 

which is larger than the peak fiux from the eddy by a factor ~ 7^ . The synchrotron fiux 
in a fixed observer band above the peak frequency is larger for the inter-eddy medium by 
an additional factor of We thus conclude that the synchrotron emission observed 

in the optical band is completely dominated by the inter-eddy medium. We therefore ignore 
eddies when we estimate the optical synchrotron fiux. 
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The situation is different for tlie IC emission. Let us write tlie synclirotron flux as seen 
by a typical electron in tlie inter-eddy medium as fgyn (this is easily estimated from the 
calculation above). The observed IC luminosity due to electrons in the inter-eddy medium 
fl^ is then 

oc arfsyuN^T^ (52) 
while the IC emission from an eddy pointing towards the observer fi^"^^ is 

f^;^''^ oc aT{ltfsyn)N,d{T^t)Vlt ~ arfsyuN^T^ (53) 

We see that the two contributions are equal. Therefore, both components in the shell fluid 
contribute equally to the gamma-ray IC flux. Of course, the inter-eddy contribution will 
vary smoothly over the duration of the burst, whereas the contribution from the eddies will 
be highly variable. 

Using these results, we may easily estimate the values of various parameters in GRB 
0803 19B corresponding to the relativistic turbulence model. Note that, since the synchrotron 
radiation, or optical flux, comes from the inter-eddy plasma, it satisfies the same equations 
as derived in §4. Moreover, the IC flux has no dependence on the Lorentz factor of turbulent 
eddies (eq. |53]). Therefore, equations fl39l) - P8|) may be directly used for the relativistic 
turbulence model provided we replace 6t with the burst duration 50 s, and take N^., Eb, 
and fic to be two times larger than the values given by these equations (the factor of two is 
to count both the eddies and the inter-eddy medium). 

Setting 6t = 50s (eq. [50]), fop = 10 Jy, = 650 keV and rji ^ = 1 (i.e., r] = Ui/ua ^ 25) 
into equation (139!) we obtain the transition radius 

Rtr = 2x 10^^ cm. (54) 

The other parameters follow from equations (H0l) - fH8l) 

r = 81i?}/^ (55) 

[254 R<Rtr, 
li=< ^ (56) 
l349i?i, R>Rtr, 

19.4x10^2^ J R<Rtr, 
N,= l (57) 

,1.3X10^3^1, R>Rtr, 
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|2.1i?;Lg , R < Rtr, 

Y=l ^ (58) 

[5.5 ni, R>Rtr, 
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(2.6 X lO^^erg) i^l6^ R < Rtr, 

Eb = { (59) 



,(1.9 X lO^^grg) R^^^ R> R 



tr ) 
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;i.3x lO^ierg) < i?^., 

= < (60) 
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,(2.6 X lO^^erg) > 

'(0.3s) Rle, R < Rtr, 

tsyn = { (61) 

(2.8s) i?S, R>R, 



37 

0.07s pi6 R ^ R 

,Lobs,52 -^16' ^ -^tr, 

io.80i?r/, i?<i?f 



(62) 



46 i -IL ^ JLtr, 

/.c-3 = <! ^ (63) 

As discussed above, we include only the inter-eddy medium for calculating the synchrotron 
component of the emission, but we include both the eddies and the inter-eddy medium 
when calculating the IC component. The synchrotron cooling time given by equation (1611) 
applies only to electrons in the inter-eddy medium; the timescale is larger by a factor of jt 
for electrons in eddies. We note that the IC flux should be increased by a factor of 2.5 to 
allow for the expansion of the source shell (Appendix A). Also, the flux will be larger by a 
factor Ned'jf/Ni than given in equation fl63|) if there are more electrons in eddies than in the 
inter-eddy medium. 



5.1. Relativistic turbulence: a consistent model for GRB 080319B 

Figure 2 is similar to Fig. 1, but shows what happens when we use the relativistic 
turbulent model. We find a very good match both with the observations and with various 
consistency conditions when we choose R ~ 10^'' cm. For this choice of R, we have (i) 
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Fig. 2. — Similar to Fig. 1, but for the relativistic turbulence model. Note that the predicted 
IC flux (solid lines in the topic left panel) is perfectly consistent with observations (horizontal 
dotted line). The shaded region indicates the range of source radii Riq that is consistent 
with constraints on the total energy (solid line in the top right panel) and the cooling time 
(solid line in the bottom left panel) constraints (see text for details). We note that the ratio 
of IC and synchrotron luminosities is larger by ~ 10 than the value of Y shown as a dashed 
line in the lower left panel (see eq. [9] for the definition of Y) due to two different effects each 
contributing a factor ~ 3; (1) the commonly used Compton-F is larger than Y in eq. [9] by 
a factor 3 due to p dependent factors, and (2) the synchrotron photon energy density inside 
the shell is larger than the naive estimate of L^yn/ (47r i?^c) by a factor ~ 3. 
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gamma-ray flux /ic-s predicted to be close to the observed flux, (ii) modest requirement 
for the total isotropic energy ~ 10^^ erg, (iii) Eb ~ -Eg/lO, i.e., approximate equipartition 
between magnetic and thermal energy, (iv) tcooi ~ 50 s, i.e., the cooling time is comparable to 
the burst duration and thus consistent with efficient radiation, (v) Y ~ fe^ft^l, i.e., consistent 
with gamma-rays dominating the emission, (vi) F ~ 250 as inferred for GRBs in general 
from a variety of observations (Lithwick & Sari, 2001), and (vii) electron Lorentz factor in 
the inter-eddy medium 7j ~ 500 and in eddies 7i/7t ~ 50, which are quite reasonable. 

The shaded band in Fig. 2 shows the range of i?i6 that is consistent with our two 
primary constraints. First, we require the total isotropic energy over the duration of the 
burst to be no larger than 10^^ erg. This constrains i?i6 to lie in the range 0.31 — 32. Second, 
we require the cooling time tcooi to lie within a factor of 10 of the burst duration 50 s. This 
gives the constraint 6.5 < Ri^ < 50. Requiring both conditions to be satisfied simultaneously 
restricts Riq to lie in the range 6.5 — 32, as shown in Fig. 2. Within this range, the predicted 
gamma-ray flux agrees remarkably well with observations. 

Note that the deceleration radius for the blast wave is 

R, = (t-^^^T^) = X 10''^^ E]l,\8tl^^sfl'n'l\ (64) 
\47™mpC(l + 2;)/ 

where n is the mean particle density of the circumstellar medium within the radius R^- 
It is interesting that i?^ lies in the middle of the allowed range for the source distance R. 
It provides independent conffimation that the prompt radiation in GRB 0803 19B was not 
produced in internal shocks - there is no reason why internal shocks should occur at the 
deceleration radius. 

In addition to the various successes described above, the relativistic turbulence model 
explains all the major qualitative features observed in the 7-ray and optical lightcurves of 
GRB 080319B during the initial ~ 10^ s, i.e., before the onset of forward shock emission. 

Since the gamma-ray emission (via IC) arises partly from eddy electrons and partly 
from the inter-eddy medium, we expect the gamma-ray lightcurve to consist of a smooth 
slowly-varying component plus a large number of sharp spikes. This is the case for most 
GRBs, including GRB 080319B (e.g.. Fig. 1 in Racusin et al. 2008). The relative fluxes 
in the two components provide information on the relative numbers of electrons in the two 



^The ratio of energies in the IC radiation and the synchrotron emission is ^ lOF. This is in part because 
the parameter y, defined in eq. [9l differs from the commonly used Compton-F parameter by a factor ^ 3 
(due to p-dependent factors not included in eq. [9]), and the mean synchrotron photon energy density in the 
shell is larger than Lsj,„/(47ri?^c) by a factor 3. 
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media. We assumed in our model (for convenience) that the numbers are roughly equal and 
this is reasonably consistent with the observations. As already mentioned, by combining the 
7-ray variability time of ~ 0.5s with the burst duration of 50s, we infer that 7t ~ 10 for 
GRB 080319B. 

Since the synchrotron emission is generated by inter-eddy electrons, the optical light 

curve is expected to be much less variable than the IC-dominated gamma-ray emission. 
This is indeed the case for GRB 0803 19B. At the same time, the overall duration of the 
optical and gamma-ray lightcurves are expected to be similar, as observed. 

The optical lightcurve of GRB 080319B showed an initial rapid rise by more than an 
order of magnitude in flux (Fig. 3 in Racusin et al. 2008), whereas the 7-ray lightcurve 
showed a much less rapid rise. This finds a natural explanation. According to the turbulent 
model of GRB 080319B, the synchrotron frequency z/j was below the optical band. Therefore, 
the optical spectrum is predicted to be very soft: F^, oc If we assume that z/j initially 

started off at a somewhat lower frequency and later settled down at a larger value, say by a 
factor of ~ 3, then the optical flux would increase by nearly a factor of 20. This explanation 
might indicate that the gamma-ray peak energy i/ic should also increase with time, whereas 
in fact decreased by a small amount (from 750 keV to 550 keV). To explain this, we would 
need to invoke that the electron Lorentz factor 7^ decreased by a factor of about 2 during 
this time. Note that the reason for the much less rapid increase of the 7-ray fiux is that the 
synchrotron peak fiux, which is proportional to the number of electrons and the magnetic 
field strength, is a slowly varying function of time. 

Another property of the rclativistic turbulence model is that we should continue to see 
emission in the 7-ray band for a time duration somewhat longer than the prompt optical 
lightcurve duration. The reason is that there is a very high probability that a few eddies 
lying a little bit outside of F^^ will point toward the observer, thereby slightly lengthening 
the burst duration in the gamma-ray band (see Fig. 1 in Narayan & Kumar, 2008). This 
effect is clearly seen in the hghtcurves of GRB 080319B; the optical LG started faUing off at 
43s whereas the steep dechne of the gamma-ray LG began at 51s. 

A prediction of the model is that the synchrotron and IG spectra should be the same. 
In particular, since our solution for GRB 080319B requires R ~ 10^'^ cm > Rtr, or i/i < 2 eV, 
the spectral index in the optical band during the burst should have been the same as the 
high energy index in the gamma-ray band, i.e. /3 = 2.87. It is unfortunate that there were no 
measurements of optical spectrum during the burst. The first measurement was at i ~ 10^ s 
when it was found that f3 = 0.55 or /j, oc z/^°'^^ (Wozniak ct al. 2008). This measurement 
would seem to call into question the prediction of the SSG model. 
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It is, however, interesting to note that the optical hghtcurve showed a sharp break at 
about 90 s. Prior to this time the flux scaled as t~^'^ and after this time the flux decreased 
as (see Fig. 2 of Kumar & Panaitescu or Racusin et al. 2008). This suggests that 90s 
marked a transition from one source of radiation to another, and that (3 = 0.55 at ~ lO^s 
corresponds to the second source which gave rise to the fi, cx part of the lightcurve and 
possibly unrelated to the prompt radiation. 

The optical lightcurve decline of t~^'^ between 43s and 90s is roughly consistent with 
the expectation of the relativistic turbulence model after the source is turned off at t ~ 43s; 
the observed radiation in this case is the large-angle emission (LAE) from photons arriving 
from angles larger than T~^, leading to a flux decline of or ~ when /5 ~ 3 (Kumar 
& Panaitescu, 2000). So the steeply declining optical lightcurve at the end of the GRB 
provides an indirect confirmation of a steep spectrum in the optical band. We note that the 
temporal behavior of flux from an adiabatically expanding shell of angular size F"^ is similar 
to the large angle emission. The decay index of the lightcurve for adiabatic expansion is 
somewhat steeper than LAE for a given f3 (Barniol-Duran & Kumar, 2008), and therefore 
this is a preferred mechanism for the observed optical flux during 43-90s. 

The gamma-ray lightcurve in the 15-150 keV at the end of the burst was seen to fall 
off even faster than the optical flux at the end of the burst. This is a puzzling behavior, 
and unlikely to be due to LAE. The reason is that the spectral index in this band was close 
to zero during the burst, and therefore the LAE flux decline should be ~ — unless the 
peak frequency fell off from 650 keV to less than 100 keV at the end of the burst which 
seems unlikely. The only natural explanation for the steep decline of the 7-ray flux is that 
the angular size of the source was ~ F~^, and gamma-rays for t > 51s were from the 
adiabatically cooling source; the IC flux from an adiabatically cooling shell declines much 
faster than the synchrotron lightcurve (Barniol-Duran & Kumar, 2008). 

What about the fall-off of the optical flux as t^^'^ for t > 10^ s? It cannot be LAE for 
the reasons described above. We offer an explanation for this part of the optical lightcurve 
that requires the GRB jet to be a Poynting outflow (see §5.2.4). A Poynting jet traveling 
outward from the center of the star cannot avoid sweeping up and accumulating some amount 
of baryonic material at its head. In subsequent jet expansion this baryonic gas is cooled, and 
at the deceleration radius it is heated once again by the reverse shock. The optical emission 
from this reverse-shock heated gas might be responsible for the lightcurve for 10^ 10^ s. 
We know from equation (168|) below that N^/Nrs ~ 10 (assuming that the kinetic energy in 
the baryonic head of the Poynting jet is of order the explosion energy). Therefore, the optical 
flux from the reverse shock is a factor ~ 10 smaller than the prompt gamma-ray source. The 
more slowly declining reverse shock flux took over from the very rapidly declining flux from 



- 25 - 



the early GRB tail at t ~ 10^ s, and continued to dominate the lightcurve until the even 
more slowly declining, but weaker, forward shock optical emission took over at t ~ 10^ s. 
We note that the effect of a narrow jet, with opening angle ~ F"^, is very weak on the 
emergent lightcurve decay for a long period of time when the jet is propagating in a medium 
with density falling off as (Kumar & Panaitescu, 2000); the density in the circumstellar 
medium of GRB 080319B is in fact inferred to be by the late time afterglow data, cf. 
Racusin et al. (2008), Kumar & Panaitescu (2008). 

5.2. Where exactly is the turbulent region located? 

As we have seen, the relativistic turbulence model gives robust estimates for various 
source parameters such as the radius, bulk Lorentz factor and energy of the shell, the number 
and typical Lorentz factors of the radiating electrons, etc. Using these results we now attempt 
to infer where the radiating region is located within the context of a dynamical model of 
GRBs. 



5.2.1. Not in internal shocks 

The internal shock model, including all reasonable variations, is firmly ruled out, as we 
have discussed in §4. Inclusion of relativistic turbulence within the context of this model 
will not salvage the situation unless we take 6t to be the burst duration (as shown in §5). 
However, in that case we are dealing with a situation in which the emission region is close 
to the deceleration radius, which is no longer an internal shock. 

The primary motivation for the internal shock model is to explain the rapid variability 
observed in the gamma-ray lightcurves of GRBs (Sari & Piran 1997). The relativistic tur- 
bulence model described in Narayan & Kumar (2008) has a completely different explanation 
for the variability. In particular, this model no longer needs to assume equation (jlj), which 
is the key relation in the internal shock model. Therefore, we see no reason to retain the 
internal shock picture. 

5.2.2. Not in the forward shock 

In the standard model of GRBs, the collision of the relativistic ejecta with the external 
medium causes a pair of shocks to be generated: a forward shock (FS) which is driven into 
the external medium and a reverse shock (RS) which is driven into the ejecta. 
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We can rule out the FS by considering the number of electrons we need for producing 
the observed radiation. From equation fl571) we see that the radiating region must have about 
6 X 10^^ electrons. However, the number Nps of electrons/protons processed in the FS must 
satisfy, by a simple energy argument, 

NFsr^mpC^ = E/2, or Nps = '^^}^!^ n p = 2 x W'^E^SM^ + ^Y' R^^ (65) 

mpC[l + z)n 

For the particular case of GRB 080319B this gives Nps ~ 10^^, which is smaller than the 
number of electrons needed by a factor ~ 10^. This is a large discrepancy, so we can discard 
the FS as the location of the relativistic turbulence. 



5.2.3. Relativistic turbulence in the reverse shock? 

Could the relativistic turbulence be located in the RS? Let the GRB ejecta be composed 
of protons and electrons, and let us take the Lorentz factor of the RS front with respect 
to the unshocked ejecta to be Trs- By applying pressure equilibrium across the contact 
discontinuity between the FS and RS fluid, we find the number of electrons Nrs that have 
been processed through the RS to be 

Nrs = Nfs^. (66) 

J- RS 

Using equations fl55l) & fl65|) . and the parameters for GRB 0803 19B, we find 

Nrs = 1.3 X 10^^E55i?i('r^^i. (67) 

Thus, the ratio of the number of electrons needed for optical/gamma-ray radiation (eq. [57j) 
and Nrs is given by 

^ = 5£-r- (68) 

If gamma-rays were to arise in the reverse shock then we expect E^,^ ~ 0.4. The reason 
is that half the energy of the blast wave is in the reverse shock at the deceleration radius, 
and this energy is efficiently radiated when the cooling frequency is close to z/j, as seems 
to be the case for GRB 080319B. Moreover, we presumably require Trs ^ 10 in order for 
the shocked gas to have a turbulent 7t ~ 10. The requirement that tcooi ~ 50 s means that 
Rn ~ 1. Therefore, we find from the above equation that N^/Nrs ~ 10. This ratio might be 
closer to unity provided that protons carry a much larger fraction of the blast wave energy, 
so that -Ess is ~ 1 — 2 rather than 0.4. 



-27- 



The interesting result that it is possible to have A^e ~ ^rs suggests that the turbulence 
is perhaps produced in the RS-heated GRB ejccta. The ratio of energies in magnetic fields 
and particle kinetic energy in this case is ~ 0.1 (fig. 2), which is similar to the value derived 
for the Crab pulsar at the wind termination shock (Kennel & Coroniti, 1984). Presumably, 
the turbulence is a natural consequence of a relativistic shock. For instance, the contact 
discontinuity surface separating the FS and RS region is known to suffer from the Rayleigh- 
Taylor instability. Could this explain the turbulence? In the shell comoving frame, the 
growth rate of the Rayleigh- Taylor instability at the interface of a relativistic RS and FS 
can be shown to be 

co^ = fkg, (69) 

where g ~ c^T/R is the effective gravitational acceleration in the shell comoving frame, 
k — 2t:IT/R is the wavenumber of the perturbation, / = 3/(8r|j5 — 5), and Trs is the 
Lorentz factor of the RS front with respect to the unshocked ejecta. For Trs » 1 the above 
equation reduces to 

u ~ ~ ^. (70) 

where 5t' is the GRB duration in the shell comoving frame. Thus, the number of c-folds 
by which the Rayleigh- Taylor mode can grow is ~ uj{5t') ~ ^^^"^^^s- '^^^ eddy scale £ of 
interest to the IC problem is R/ (FF^^^) or £ ~ F^j^. Perturbations on this scale will undergo 
F^y^ e-folds of growth, i.e., the amphtude increases by less than a factor 2. Therefore, 
the Rayleigh-Taylor instability is not sufficiently potent to generate the highly relativistic 
turbulence we need. 

Recently Goodman & MacFadyen (2007) and Milosavljevic, Nakar & Zhang (2007) have 
discovered interesting instabilities, resulting from a clumpy circumstellar medium and an 
initially anisotropic blastwave respectively, which lead to vorticity generation downstream of 
the shock front. These instabilities have been further studied by Sironi & Goodman (2007), 
and Milosavljevic et al. (2007) to investigate the generation of magnetic fields in relativistic 
shocks. Couch, Milosavljevic & Nakar (2008) have found another instability that generates 
vorticity down stream of a shockfront even when the circumstellar medium is homogeneous 
and the blastwave isotropic. We have estimated the growth rate of these instabilities and 
find that these too fail to give rise to relativistic turbulence. 

Of course, we cannot rule out the possibility that there might be other as yet unknown 
instabilities that might give rise to relativistic turbulence. Therefore, we are unable to 
discard the possibility that the prompt GRB emission originates in the RS. 
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5.2.4- Relativistic turbulence in the Poynting- dominated jet 

A Poynting-dominated jet would have a weak reverse shock (Kennel & Coroniti, 1984; 
Zhang & Kobayashi, 2005) and would not be consistent with the proposal considered in 
the previous subsection. On the other hand, such a jet probably undergoes various plasma 
instabilities at the deceleration radius. These instabilities would stir up the fluid into a state 
consistent with our model of relativistic turbulence. The instabilities would presumably heat 
up the electrons until quasi-equipartition is achieved, consistent with the results shown in 
Fig. 2. 

According to equations (15^ & (p7]l . in our model Es/Ee ~ 0.1 at -R ~ 10^^ cm. However, 
this does not rule out the Poynting outflow model. The reason is that in all of our formulae 
B really stands for the projection of the magnetic field vector perpendicular to the electron 
momentum vector i.e., B sin a where a is the pitch angle between the electron momentum and 
the magnetic field direction. For a random distribution of particle pitch angle the difference 
between B and B sin a is order unity. However, when electrons have a non-zero average 
momentum along the local magnetic field (as might be the case for particles accelerated in 
reconnection regions), the difference can be large. For instance, when the average a is 0.3 
the energy in magnetic fields is larger than that in equation (159|) by a factor ~ 10, making 
the model consistent with equipartition. 

Lyutikov and Blandford (2003) have suggested that the dissipation of magnetic energy in 
a Poynting flux dominated jet should occur at a distance of ~ 3x lO^^cm due to current driven 
instabilities (see Lyutikov 2006 for a concise summary of the model, and for a comparison 
with the baryonic outflow model). Acceleration of electrons (and positrons), and plasma 
bulk flow along the magnetic field lines at roughly the local Alfven speed are expected in the 
process of magnetic field decay/reconnection. These expectations of the Poynting outfiow 
model are roughly consistent with our findings for GRB 080319B: emission generated at 
R ^ lO^^cm and turbulent velocity field with Lorentz factor 7^ ~ 10. However, the reason 
for a very soft particle spectrum, p ~ 5, is unclear (at least to us); numerical simulations 
of particle acceleration in reconnection regions generally find a hard particle spectrum (e.g., 
Larrabee et al. 2003). 

Moreover, it is also not clear how jt and 7, should be related to the bulk P of the 
pre-instability jet. Nor is it clear why the typical Lorentz factor of electrons should be a 
modest value 7j ~ 500 with the kind of powerful accelerator one might expect in magnetic 
reconnections (other than the fact that it is energetically impossible to accelerate a large 
number of electrons, of order 10^^, to an average Loretnz factor much larger than ~ 500). 
Further investigation is required to address these questions. 
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6. Summary 

We have shown in this paper that the gamma-ray and optical data for GRB 0803 19B 
rule out the popular internal shock model for generation of the prompt radiation. According 
to this model, the duration (< 1 s) of spikes in the gamma-ray lightcurve sets an upper bound 
on the quantity i?/(2cr^), where R is the radius of the source relative to the center of the 
explosion and F is the bulk Lorentz factor. When we apply this condition, we find that it is 
impossible to fit the observed optical and gamma-ray flux simultaneously. Specifically, any 
model that fits the optical flux under predicts the gamma-ray flux by nearly two orders of 
magnitude (Fig. 1). This is an unacceptably large discrepancy which cannot be eliminated 
with any reasonable modification of the internal shock model. 

An equally powerful qualitative argument against the internal shock model is the fact 
that we find the radius R of the source to be constrained quite tightly by the observations. 
The energy required in magnetic field increases very rapidly as we decrease R, Eb oc 
whereas the energy in particles increases rapidly, E^, oc i?^^/^. Also, the coolin g t ime of 
electrons becomes too short to be compatible with observations if i? < 3 x 10^^ cmo All of 
these factors together constrain the location where the prompt emission in GRB 0803 19B 
was produced to lie within a narrow range of radius: 4 x 10^^ < i? < 8 x 10^^ cm (see Fig. 1). 
This is problematic for the internal shock model. According to this model, there is a large 
number of internal shocks among independent ejecta, with a separate shock producing each 
of the ~ 50 spikes in the gamma-ray lightcurve of GRB 080319B. Why would all the ejecta 
collide within such a narrow range of radius? Moreover, why should the radius be so close to 
the deceleration radius Rd ~ 10^'' cm, where the ejecta meet the external medium and begin 
to slow down? This coincidence is suspicious. 

All of these problems are eliminated if we give up the internal shock model and consider 
instead a model in which the variability in the gamma-ray lightcurve is produced by rela- 
tivistic turbulence in the source with random eddy Lorentz factors 7f ~ 10. In this model, 
the quantity R/{2cT'^) is no longer constrained to be less than 1 s, but only needs to be 
comparable to the burst duration ~ 50 s (Narayan & Kumar 2008). With this modification, 
we find that we obtain a remarkably consistent model of GRB 0803 18B (see Fig. 2) in which 
the prompt optical emission was produced by synchrotron emission and the gamma-rays 



^^CoUisions at a smaller radius would produce a weak optical flash with flux decreasing roughly as R^^/^^. 
The electrons would undergo very rapid cooling and produce a low energy spectrum in the gamma-ray band 
of fv cx Prompt optical observations of GRB 080319B show variations in the optical flux by less than 

a factor two for much of the 50 s duration of the burst except at the beginning and the end. Moreover, the 
low energy spectral index for the gamma-ray emission was greater than throughout the burst. 
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were the result of inverse Compton scattering. The predicted gamma-ray flux is perfectly 
compatible with observations. Also, estimates of various quantities such as the total energy, 
cooling time, Lorentz factor, etc. are all very reasonable and consistent (§5.1). The radius of 
the source is calculated to be in the range 6 x 10^^ < i? < 3 x 10^'' cm; if we select a nominal 
value R ~ 10^^ cm, we obtain an excellent fit to all the observations. 

In the context of a physical model, the picture that emerges from this model is that the 
energy of the relativistic jet in GRB 080319B was converted to optical & 7-ray radiation 
either via a relativistic reverse shock when ejecta (composed of p'^s and e~s) ran into the 
circumstellar medium or that much of the jet energy was in magnetic field that was dissipated 
close to the deceleration radius. Theoretically, it is difficult to understand how a reverse shock 
might produce relativistic turbulence with 7^ ~ 10 (§5.2.3). Also, it is easier to understand 
the optical data for the time period 10^ < t < 10'^ s if we assume a Poynting jet (§5.2.4). For 
these reasons we have a mild preference for the Poynting-dominated jet model. 

A potential problem for the Poynting jet model is that the ratio of magnetic to particle 
kinetic energy is about 0.1 for our best solution (Fig. 2). However, this ratio is similar to 
that inferred for the pulsar wind termination shock for the Crab pulsar (Kennel & Coroniti, 
1984). Moreover, this ratio of 0.1 does not rule out the Poynting model for another reason 
which is that, in all of our formulae, B is the projection of the magnetic field perpendicular 
to the electron momentum vector. Thus, if electron momenta are preferentially parallel to 
the magnetic field, then the true Eb would be larger than our estimate (easily by a factor 
10 compared to the value given in eq. and we can have Es/Ef. ~ 1. Note that electrons 
are accelerated parallel to the magnetic field in reconnection regions and so this possibility 
is not as arbitrary as it might appear. 

In the relativistic turbulence model, fiuctuations in the observed gamma-ray lightcurve 
are produced as a result of random relativistic variations in the velocity field of the source, 
with turbulent Lorentz factor 7^ ~ 10. The model predicts that there should be ~ ~ 100 
spikes in the gamma-ray light curve (Narayan & Kumar 2008), which is consistent with the 
~ 50 spikes seen in GRB 080319B. The optical synchrotron fiux is dominated by the inter- 
eddy medium rather than eddies. Therefore, we expect much less variability in the optical 
fiux, as was indeed observed. 

The model can explain the sharp rise in the optical fiux of GRB 080319B at the beginning 
of the burst. For this, we must postulate that the synchrotron peak frequency increased from 
~ 0.5 eV to ~ 1.5 eV during the first ~ 15s. Since the synchrotron peak frequency z/j was 
below the optical band (2eV), the optical fiux was smaller than the peak synchrotron fiux by 
a factor (4/z/j)^'* (the spectral index above the peak is known from gamma-ray observations). 
A modest increase in z/, by a factor of 3 early in the burst would thus produce a factor of 
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20 increase in the observed optical flux. The reason that the gamma-ray flux increased by a 
much smaller factor during the same time is that the peak IC flux is proportional to TeN^B, 
which would change little during this time period. 

The end of the gamma-ray prompt emission phase occurred ~ 8 s after the prompt op- 
tical in GRB 080319B. This is probably a result of inverse Compton emission from turbulent 
eddies lying a bit outside of the primary 1 /F cone, but which happened to point toward us 
because of a fortuitous alignment of their turbulent velocity (the probability for this hap- 
pening is of order unity). Since much of the synchrotron emission comes from non-turbulent 
fluid in between eddies, the optical flux would not have a similar effect. 
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A. Possible errors in the calculation of IC flux 



We discuss in this appendix possible sources of error in our calculation of the IC flux, 
i.e. errors associated with eq. (HHj) : according to this equation the theoretically calculated 
gamma-ray flux is smaller than the observed value by a factor ~ 20. 

A possible source of error might arise from our assumption of a homogeneous source, and 
we need to estimate its effect on the IC flux. The synchrotron peak flux is a linear function 
of magnetic field strength and the total number of electrons in the source, and therefore 
clumping of electrons and B, to lowest order, have little effect on the emergent flux. The IC 
flux is, however, affected by clumping of electrons, and we estimate the magnitude of this 
effect. 

Let us consider an extreme form of inhomogeneity where all electrons are concentrated 
in Mc clumps of each size r^. The number density of electrons in the clumps is n^, and 
the density averaged over the source volume, ~ i?^, is ng; there are no electrons in between 
clumps. Let us assume that the synchrotron power from each electron is p^. In this case the 
synchrotron luminosity of the source is 



and is independent of electron clumping. The IC luminosity depends on synchrotron flux in 
the vicinity of electrons (in clumps). The synchrotron flux is 



(Al) 




1 + 



TIqR 



(A2) 



i?2 



The IC luminosity is obtained from the above flux: 



lM^^aTLsyn{y)noR 1 + 



UqR 



(A3) 



This preprint was prepared with the AAS IM^TpC macros v5.2. 
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Or 



(A4) 



where /c = ricrlMc/ (no-R^) is the fraction of the shell volume occupied by clumps. We 
see from the above equation that IC flux can be enhanced by clumping of electrons. For 
instance, consider an example where fc = 0.1 and Mc = 1 (all electrons are in a single small 
clump). The IC flux in this case is a factor 4.5 larger than when electrons are uniformly 
distributed. This flux enhancement is about an order of magnitude smaller than what is 
needed to explain the observed flux for GRB 080319B (eq. I48ll . An even more extreme 
case of clumping could bridge the gap, however the efficiency for converting jet energy to 
radiation is very small when -C 1 as pointed out by Sari and Piran (1997). Furthermore, 
another serious problem is that a high degree of clumping leads to an increase of Ua (as 
discussed below), and that makes the SSC spectrum below the peak inconsistent with the 
observed data for GRB 080319B - unless we place the source at a distance from the center 
of explosion that is larger than the deceleration radius. 

The dependence of fic on r] = Vi/va is fairly strong and so we need to discuss the 
uncertainty in rj. We have taken rj = 25 (log 77 = 1.4), which is guided by the Konus-Wind 
low energy spectrum of f^, oc z/°-^ in the energy band 20-650 keV. This low energy spectral 
index suggests that Ua (the self-absorption frequency) should be <20 keV, and thus ?7>32; 
therefore, 77 = 25 is a conservative choice for GRB 080319B. However, is it possible that 
Ua has been overestimated in our calculation by our assumption of a homogeneous source? 
If Ua were to be smaller by a factor ~ 6 than given by equation (ITBi) then that would lead 
to a larger IC flux by factor 30 (see eq. US]), and thereby reconcile the observed and the 
theoretically expected gamma-ray flux. We show that inhomogeneities in the source cannot 
decrease z/q as long as the optical flux we observed during the burst is produced in the source. 

We calculate synchrotron self-absorption frequency {ua) when B, 7j & uq are allowed 
to vary, arbitrarily, across the source; the electron distribution is taken to be dn/dj = 
^o{l/li)~^ for 7 > 7i. Spatial variations in B, 7^ & no are subject to constraints that the 
optical flux and the IC peak frequency should be equal to the observed values. 

Our starting point is equation (6.52) of Rybicki & Lightman (1979) for the synchrotron 
absorption coefficient, a^,; f dr' is the optical depth for absorbing synchrotron photons 
of frequency v. We can show that for a power-law electron distribution and for v < Vi (the 
case of interest for 080319B) 




(A5) 
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where 



47rmeC 



6 is the angle between magnetic field and electron velocity vector, and prime denotes fre- 
quency in the source comoving frame. The synchrotron self- absorption frequency is deter- 
mined from the equation 



/ 



3^/2(p + 2)(p- l)g3 f noBsm6 



dr'^K - :.j:jrZ.J::2 l dr' '-^^^^^ ( ^ ) =1. (AT) 



Or 

,5/3^, 3V2(p + 2)(p-l)g3 f ^^^noBsm5 



1/3 



The z/^ given by the above equation is self- absorption frequency along one line of sight. Since 
an observer receives photons from an area ~ 7ri?^/r^, we should average Ua over this area. 
This average frequency is given by 

,5/3. 3V2(p + 2)(p-l)g3 r2 f^^^,noBsm5 



^ ^ 167r2m,(p + 2/3) vr/?^ j ^ ^ ^/i/a^^" ^^^^ 

Since the optical flux [fop) is proportional to / d^id tiqB sm.5 / p'^^^ , when z/^' lies above 
the optical band, it is convenient to define a new variable x = n^B sm.5 /p'^^^ , and rewrite 
the equation for v'^ using this new variable 

K'^')oc-^|rfV|. (AlO) 

The minimum of v'^^ can be obtained by requiring that 5v'^ = for an infinitesimal variation 
of x(x'), i.e. 



subject to the condition that 



d^x.' ; = 0, (All) 

y"rfV5x(x') = 0. (A12) 

The variational integral implicitly assumes that we are solving for 7j(x') which can be an 
arbitrary function as long as the IC spectral peak of the radiation emergent from the source 
matches the observed value, i.e. 1 <^ 7i ^ oo. Let us consider a special form of 5x(x') 
that is nonzero in two spherical regions of infinitesimal radius centered at x'-^ and x'2. It is 
required that 5x(x'i) = — (5x(x'2) in order to satisfy the optical flux constraint. Substituting 
this into equation (lAlip leads to 7j(x'i) = 7j(x'2). Since, x'l and x'2 were arbitrary points 
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in the source, we conclude that 7^ does not vary across the source when z/^ is minimized. 
Therefore, we can take 7^ outside of the integral in equation flAlOp . and find see that the 
minimum value of {va) is fixed by the observed optical flux. In other words, the assumption 
of a homogeneous source used in our calculations in §3 & §4 gives the smallest possible value 
for 7] or the largest IC flux. For a high degree of clumping of electrons in the source, we 
considered above, Va is larger when our line of sight passes through a clump and especially 
when only one clump lies on our sight line; there is, however, little change to as we have 
shown above. In this case an observer will receive radiation from one clump at a time, and 
will find Va larger than the case of a homogeneous shell. A larger Va for a clumpy source goes 
in the opposite direction to what we need to increase the IC flux, and this largely reverses 
the gain to the IC flux found above (eq. A4). 

We assumed in our derivation that v[ is above the optical band. A similar proof for 
the minimum of v'^ can be carried out when v[ is below the optical frequency. Moreover, we 
ignored variations of F across the source. This approximation is justified since large variation 
in F are smoothed out in less than one dynamical time. 

There is one effect that we have not included in our calculation which lowers the value 
of the self- absorption frequency Va somewhat. As synchrotron photons propagate outward 
they move though a medium where the electron density is decreasing with time (due to the 
outward expansion of the shell). As a result, a calculation based on a stationary source 
overestimates the optical depth to Thomson scattering and i>a by factors of ~ 2 & ~ 2'^/^, 
respectively. This can be seen by considering a homogeneous shell with electron density no 
located at a distance Rq from the center of explosion, with a radial thickness Rq/T'^. (The 
comoving frame shell thickness of Ro/T is obtained by causality considerations, since Rq/cT 
is the time elapsed in the shell comoving frame.) As seen by a lab frame observer, a photon 
moving outward in the radial direction takes a time 2Rq/c to cross the shell, and during this 
time the shell has moved to a larger radius and the density has decreased. For a hot shell, 
the radial width too increases with time, and the photon transit time is a bit larger still. 
A straightforward calculation of the optical depth to Thomson scattering when the density 
changes as l/R"^ shows that the shell optical depth is about half of what it is for a stationary 
medium of the same thickness and density. Substituting this optical depth into equation 
flASp for Ua we see that the synchrotron-self absorption frequency is smaller for an expanding 
shell by a factor of ~ 2^/^. 

The net effect is that the parameter 77 which is defined in equation (ITTI) should be 
reduced by a factor of 2^/^. In the case of GRB 0803 19B, we estimated ~ 25 from the 
observations. We should use a smaller value 77 ~ 16 in our formulae to obtain more accurate 
numerical estimates of quantities. This causes the predicted gamma-ray flux to be increased 
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by a factor of 2.5. 



